Population clustering as a signal for de-confinement 
in accreting compact stars 
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ABSTRACT 

We study the evolution of the rotation frequency for accreting compact stars. The discon- 
tinuous change of the moment of inertia of a rapidly rotating star due to the possible quark 
core appearance entails a characteristic change in the spin evolution. Numerical solutions have 
been performed using a model equation of state describing the deconfinement phase transition. 
Trajectories of spin evolution are discussed in the angular velocity - baryon number plane (phase 
diagram) for different accretion scenarios defined by the initial values of mass and magnetic 
field of the star, as well as mass accretion rate and magnetic field decay time. We observe a 
characteristic increase in the waiting time when a configuration enters the quark core regime. 
Overclustering of the population of Z sources of LMXBs in the phase diagram is suggested as a 
direct measurement of the critical line for the deconfinement phase transition since it is related 
to the behaviour of the moment of inertia of the compact star. 



Subject headings: accetion, accretion disks - 
— pulsars: general — X-rays: binaries 

Introduction 



stars: evolution — stars: interiors — stars: magnetic fields 



The deconfinement phase transition into a 
plasma of quasifree quarks and gluons is expected 
to occure under conditions of sufficiently high tem- 
peratures and/or densities as, e.g., in heavy-ion 
collisions, a few microseconds after the big bang 
or in the cores of superdense stellar objects. 

For the latter, signals of a phase transition 
have been suggested in the form of characteris- 
tic changes of observables (Weber 1999), such 
as the pulse timing (Glendenning, Pei & Weber 
1997; Chubarian et al. 2000), brightness (Dar & 
DeRiijula 2000) and surface temperature (Schaab 
et al. 1996; Blaschke, Grigorian & Voskresensky 
2000) of the isolated pulsars during their evolu- 
tion. 
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Recently, quasiperiodic brightness oscillations 
(QPOs) in low-mass X-ray binaries (LMXBs) have 
been observed (van der Klis 2000) which provide 
new mass and radius constraints for compact ob- 
jects (Miller, Lamb & Psaltis 1998). According 
to the model of Stella and Vietri (Stella & Vi- 
ctri 1999; Stella, Vietri & Morsink 1999) these 
constraints would allow only very compact objects 
made up of strange quark matter (Li et al. 1999). 
Due to the mass accretion flow these systems are 
candidates for the most massive compact stars at 
the limit to the formation of black holes. There- 
fore, if the possible deconfinement phase transition 
in compact stars exists at all, we expect it to oc- 
cur in accreting LMXBs (Chubarian et al. 2000; 
Glendenning & Weber 2000). 

In this paper we study the spin evolution of 
compact stars with possible quark matter cores 
for a disc accretion model and show that a pro- 
nounced signal of a deconfinement transition in 
the star interior can be obtained for accretors 
with long-lived magnetic fields and large angular 
momentum transfer. The so called Z sources of 
QPOs in LMXBs are suggested as objects which 
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should predominantly populate the region of the 
suspected phase border between hadronic stars 
and quark core stars (QCSs) (Chubarian et al. 
2000). We suggest this population clustering as 
a signal for deconfined quark matter in accreting 
compact stars. 

2. Evolution in the — N diagram 

Our investigation is based on the classifica- 
tion of rotating compact star configurations in the 
plane of angular velocity f2 and baryon number TV, 
the so called phase diagram (Blaschke, Grigorian 
& Poghosyan 2000). 

We consider in our calculations a generic form 
of a relativistic equation of state with a deconfinc- 
ment transition (Chubarian et al. 2000), which 
describes the hadronic matter by a Walecka model 
and the quark matter by a dynamical confinig 
model (Blaschke et al. 1999; Blaschke & Tandy 
2000). The deconfinement phase transition is ob- 
tained by imposing Gibbs' conditions for phase 
equilibrium with the constraints that baryon num- 
ber as well as electric charge of the system are 
globally conserved (Glendenning 1992). The on- 
set of deconfinement (mixed phase) is obtained at 
a density of about 1.5 times the nuclear saturation 
density. 

It has been shown that the critical line N cr a (SI) , 
which separates the region of QCSs from the 
hadronic ones, is correlated with the local maxima 
of the moment of inertia with respect to changes 
of the baryon number at given SI due to the change 
of the internal structure of the compact object at 
the deconfinement phase transition. Therefore, we 
expect that the rotational behavior of these ob- 
jects changes in a characteristic way when this 
line is crossed. The consequence will be an in- 
crease of the population of stars at this critical 
line which could be observed, provided that a suf- 
ficiently large number of accretors will be detected 
with their masses and spin frequencies (Blaschke, 
Grigorian & Poghosyan 2000). 

We consider the spin evolution of a compact 
star under mass accretion from a low-mass com- 
panion star as a sequence of stationary states 
of configurations (points) in the phase diagram 
spanned by SI and N. The process is governed 
by the change in angular momentum of the star 



-(I(N,Q) Sl)=K cxt , 
at 



where 



K eilt = dGMM 2 r -N out 



(1) 



(2) 



is the external torque due to both the specific 
angular momentum transfered by the accreting 
plasma and the magnetic plus viscous stress given 
by N out = nn 2 r- 3 , k = 1/3 (Lipunov 1992). For 
a star with radius R and magnetic field strength 
B, the magnetic moment is given by /x = R 3 B. 

1/3 

The co-rotating radius r c = (GM/S1 2 ) is very 
large (r c 3> r ) for slow rotators. The inner radius 
of the accretion disc is 
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where \i c is that value of the magnetic moment 
of the star for which the disc would touch the 
star surface. The characteristic Alfvcn radius for 
spherical accretion with the rate M = mN is 

rA = i2fi~ 4 GMM 2 ) . Since we are interested 
in the case of fast rotation for which the spin-up 
torque due to the accreting plasma in Eq. (2) is 
partly compensated by N out , eventually leading to 
a saturation of the spin-up, we neglect the spin- 
up torque in N out which can be important only for 
slow rotators (Ghosh & Lamb 1979), 

From Eqs. (1), (2) one can obtain the first order 
differential equation for the evolution of angular 
velocity 



dn = gextC/V, tl)-K int (N, SI) 

dt i(N,n) + n(di(N,n)/dfi) N ' 



(3) 



where 



K int {N, Si) = SlN(dI(N, Sl)/dN) a . (4) 

Solutions of (3) are trajectories in the SI — N 
plane describing the spin evolution of accreting 
compact stars, see Fig. 1. Since I(N,Si) exhibits 
characteristic functional dependences (Blaschke, 
Grigorian & Poghosyan 2000) at the deconfine- 
ment phase transition line N CI - lt (Sl) we expect ob- 
servable consequences in the P — P plane when 
this line is crossed. 
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In our model calculations we assume that 
both the mass accretion and the angular mo- 
mentum transfer processes are slow enough to 
justify the assumption of quasistationary rigid 
rotation without convection. The moment of 
inertia of the rotating star can be defined as 
I(N,Q.) = J(N,Q)/Q , where J(N,Q) is the an- 
gular momentum of the star. For a more detailed 
description of the method and analytic results we 
refer to (Chubarian et al. 2000) and the works of 
Hartle (1967); Hartle & Thorne (1968), as well 
as Sedrakian & Chubarian (1968a, b). 

The time dependence of the baryon number for 
the constant accreting rate N is given by 

N{t) = N(t ) + (t - t )N . (5) 

For the magnetic field of the accretors we consider 
the exponential decay (Bhattacharya & van den 
Hcuvel 1979) 

B{t) = [B(0) - Boo] exp(-i/r B ) + B x . (6) 

We solve the equation for the spin-up evolution 
(3) of the accreting star for decay times 10 7 < 
Tslyr] < 10 9 and initial magnetic fields in the 
range 0.2 < 5(0) [TG] < 4.0. The remnant mag- 
netic field is chosen to be B^ = 10 _4 TG 1 (Page, 
Geppert & Zannias 2000). 

At high rotation frequency, both the angular 
momentum transfer from accreting matter and the 
influence of magnetic fields can be small enough to 
let the evolution of angular velocity be determined 
by the dependence of the moment of inertia on 
the baryon number, i.e. on the total mass. This 
case is similar to the one with negligible magnetic 
field considered in (Shapiro & Tcukolsky 1983; 
Burderi et al. 1999; Chubarian et al. 2000) 
where fj, < /j, c in Eq. (3), so that only the so 
called internal torque term (4) remains. 

In Fig. 1 we show evolutionary tracks of ac- 
cretors in phase diagrams (left panels) and show 
the corresponding spin evolution Q(t) (right pan- 
els). In the lower panels, the paths for possi- 
ble spin-up evolution are shown for accretor mod- 
els initially having a quark matter core (N(0) = 
1.55 Nq, O(0) = 1 Hz). The upper panels show 
evolution of a hybrid star without a quark matter 
core in the initial state (N(0) = 1.4 N Q , Q(0) = 1 

H TG= 10 12 G 



Hz), containing quarks only in mixed phase. We 
assume a value of N corresponding to observa- 
tions made on LMXBs, which are divided into Z 
sources with N <~ 10^ 8 A r Q /yr and A(toll) sources 
accreting at rates N <~ lO _lo iV0/yr (Glendcnning 
& Weber 2000; Bhattacharya & van den Heuvel 
1979; van der Klis 2000). 

For the case of a small magnetic field decay time 
tb = 10 7 yr (solid and dotted lines in Fig.l) the 
spin-up evolution of the star cannot be stopped by 
the magnetic braking term so that the maximal 
frequency consistent $l ma x (N) with stationary ro- 
tation can be reached regardless whether the star 
did initially have a pure quark matter core or not. 

For long lived magnetic fields (tb = 10 9 yr, 
dashed and dot-dashed lines in Fig.l) the spin-up 
evolution deviates from the monotonous behaviour 
of the previous case and shows a tendency to sat- 
urate. At a high accretion rate (dot-dashed lines) 
the mass load onto the star can be sufficient to 
transform it to a black hole before the maximum 
frequency could be reached whereas at low accre- 
tion (dashed lines) the star spins up to the Kepler 
frequency limit. 

3. Waiting time and population clustering 

The question arises whether there is any charac- 
teristic feature in the spin evolution which distin- 
guishes trajectories that traverse the critical phase 
transition line from those remaining within the ini- 
tial phase. 

For an accretion rate as high as N = 10~ 8 Nq/ 
yr the evolution of the spin frequency in Fig. 1 
shows a plateau where the angular velocity re- 
mains within the narrow region between 2.1 < 
f2[kHz] < 2.3 for the decay time tb = 10 9 yr 
and between 0.4 < 0[kHz] < 0.5 when t b = 10 7 
yr. This plateau occurs for stars evolving into 
the QCS region (upper panels) as well as for stars 
remaining within the QCS region (lower panels). 
This saturation of spin frequencies is mainly re- 
lated to the compensation of spin-up and spin- 
down torques at a level determined by the strength 
of the magnetic field. In order to perform a more 
quantitative discussion of possible signals of the 
deconfinement phase transition we present in Fig. 
2 trajectories of the spin-up evolution in the P — P 
plane for stars with N(0) = 1.4 N Q and fi(0) = 1 
Hz in the initial state; the four sets of accretion 
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rates and magnetic field decay times coincide with 
those in Fig. 1. 

When we compare the results for the above 
hybrid star model (solid lines) with those of a 
hadronic star model (quark matter part of the 
hybrid model omitted; dotted lines) we observe 
that only in the case of high accretion rate (N = 
10~ 8 -/V Q /yr, e.g. for Z sources) and long-lived 
magnetic field (tb = 10 9 yr) there is a significant 
difference in the behaviour of the period deriva- 
tives. The evolution of a star with deconfinement 
phase transition shows a dip in the period deriva- 
tive in a narrow region of spin periods. This fea- 
ture corresponds to a plateau in the spin evolu- 
tion which can be quantified by the waiting time 
t = P/P = fi/fi. In Fig. 3 (lower and mid- 
dle panels) we present this waiting time in de- 
pendence on the rotation frequencies v = fl/(2n) 
for the relevant case labeled (9,-8) in Figs. 1,2. 
The comparison of the trajectory for a hybrid star 
surviving the phase trasition during the evolution 
(solid line) with those of a star evolving within the 
hadronic and the QCS domains (dotted line and 
dashed lines, respectfuly) , demonstrates that an 
enhancement of the waiting time in a narrow re- 
gion of frequencies is a characteristic indicator for 
a deconfinement transition in the accreting com- 
pact star. 

The position of this peak in the waiting time 
depends on the initial value of magnetic field, see 
Fig. 3. In the middle and lower panels of that 
Figure, we show the waiting time distribution for 
B(0) = 0.75 TG and B(0) = 0.82 TG, respec- 
tively. Maxima of the waiting time in a certain 
frequency region have the consequence that the 
probability to observe objects there is increased 
(population clustering). In the upper panel of 
this Figure the spin frequencies for observed Z 
sources in LMXBs with QPOs (van der Klis 2000) 
are shown for comparison. In order to interprete 
the clustering of objects in the frequency inter- 
val 225 < v [Hz] < 375 as a phenomenon related 
to the increase in the waiting time, we have to 
chose initial magnetic field values in the range 
1.0 > B(0)[TG] > 0.6 for the scenario labeled (9,- 
8), see also the dashed lines in Fig. 2. 



4. Signal for Deconfinement in LMXBs 

The results of the previous section show that 
the waiting time for accreting stars along their 
evolution trajectory is larger in a hadronic con- 
figuration than in a QCS, after a time scale when 
the mass load onto the star becomes significant. 
This suggests that if a hadronic star enters the 
QCS region, its spin evolution gets enhanced thus 
depopulating the higher frequency branch of its 
trajectory in the fl — N plane. 

In Fig. 4 we show contours of waiting time 
regions in the phase diagram. The initial baryon 
number is N(0) — 1.4A Q and the initial magnetic 
field is taken from the interval 0.2 < B(0)[TG] < 
4.0 . 

The region of longest waiting times is located in 
a narrow branch around the phase transition bor- 
der and does not depend on the evolution scenario 
after the passage of the border, when the depop- 
ulation occurs and the probability to find an ac- 
creting compact star is reduced. Another smaller 
increase of the waiting time and thus a popula- 
tion clustering could occur in a region where the 
accretor is already a QCS. For an estimate of a 
population statistics we show the region of evolu- 
tionary tracks when the values of initial magnetic 
field are within 0.6 < B(0)[TG] < 1.0 as suggested 
by the observation of frequency clustering in the 
narrow interval 375 > z/[Hz] > 225, see Fig. 3. 

As a strategy of search for QCSs we suggest to 
select from the LMXBs exhibiting the QPO phe- 
nomenon those accreting close to the Eddington 
limit (Bhattacharya & van den Heuvel 1979) and 
to determine simultaneously the spin frequency 
and the mass (Lamb & Miller 2000) for sufficiently 
many of these objects. The emerging statistics of 
accreting compact stars should then exhibit the 
population clustering shown in Fig. 4 when a de- 
confinement transition is possible. If a structure- 
less distribution of objects in the — N plane will 
be observed, then no firm conclusion about quark 
core formation in compact stars can be made. 

For the model equation of state on which the 
results of our present work are based, we expect a 
baryon number clustering rather than a frequency 
clustering to be a signal of the deconfinement tran- 
sition in the compact stars of LMXBs. The model 
independent result of our study is that a pop- 
ulation clustering in the phase diagram for ac- 
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creting compact stars shall measure the critical 
line N CT i t (Cl) which separates hadronic stars from 
QCSs where the shape of this curve can discrimi- 
nate between different models of the nuclear EoS 
at high densities. 
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Fig. 1. — Spin evolution of an accreting compact 
star for different decay times of the magnetic field 
and different accretion rates. Upper panels: ini- 
tial configuration with N(0) — 1.4 Nq; Lower 
panels: N(0) = 1.55 iV Q ; Q(0) = 1 Hz in both 
cases. The numbers in the legend box stand for 
(log(r B [yr]), log(iV[JV /yrJ). For instance (9,-8) 
denotes t b = 10 9 yr and N = 10" 8 7V Q / yr. 
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Fig. 2. — Spin-up evolution of accreting compact 
stars in the P — P diagram. For labels and ini- 
tial values see to Fig. 1. The region of the ver- 
tical dashed lines corresponds to the clustering 
of periods observed for LMXBs with QPOs. A 
dip (waiting point) occurs at the dcconfincmcnt 
transition for parameters which correspond to Z 
sources (N = 10 _8 A r Q /yr) with slow magnetic 
field decay (tb — 10 9 yr). 



Fig. 3. — Upper panel: frequency interval for ob- 
served Z source LMXBs (van der Klis 2000); mid- 
dle panel: waiting times r = P/P for scenario 
(9,-8) and initial magnetic field B(0) = 0.75 TG; 
lower panel: same as middle panel for -B(O) = 0.82 
TG. Spin evolution of a hybrid stars (solid lines) 
shows a peak in the waiting time characteristic for 
the deconfinement transition. Hadronic stars (dot- 
ted lines) and QCSs (dashed lines) have no such 
structure. 
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Fig. 4. — Regions of waiting times in the phase di- 
agram for compact hybrid stars for the (9,-8) sce- 
nario. For an estimate of a population statistics 
we show the region of evolutionary tracks when the 
interval of initial magnetic field values is restricted 
to 0.6 < B(0)[TG] < f.0. Note that the probabil- 
ity of finding a compact star in the phase diagram 
is enhanced in the vicinity of the critical line for 
the deconfinement phase transition 7V cr i t (r2) by at 
least a factor of two relative to all other regions in 
the phase diagram. 
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